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Abstract — This paper considers the problem of robust stabil- 
ity for a class of uncertain quantum systems subject to unknown 
perturbations in the system Hamiltonian. Some general stability 
results are given for different classes of perturbations to the 
system Hamiltonian. Then, the special case of a nominal linear 
quantum system is considered with either quadratic or non- 
quadratic perturbations to the system Hamiltonian. In this case, 
robust stability conditions are given in terms of strict bounded 
real conditions. 

I. Introduction 

An important concept in modern control theory is the 
notion of robust or absolute stability for uncertain nonlinear 
systems in the form of a Lur'e system with an uncertain 
nonlinear block which satisfies a sector bound condition; 
e.g., see [1]. This enables a frequency domain condition for 
robust stability to be given. This characterization of robust 
stability enables robust feedback controller synthesis to be 
carried out using H°° control theory; e.g., see [2]. The aim 
of this paper is to extend classical results on robust stability 
to the case of quantum systems. This is motivated by a desire 
to apply quantum H°° control such as presented in [3], [4] 
to nonlinear and uncertain quantum systems. 

In recent years, a number of papers have considered the 
feedback control of systems whose dynamics are governed 
by the laws of quantum mechanics rather than classical 
mechanics; e.g., see [3]— [15]. In particular, the papers [12], 
[16] consider a framework of quantum systems defined in 
terms of a triple (S,L,H) where S is a scattering matrix, 
L is a vector of coupling operators and H is a Hamiltonian 
operator. The paper [16] then introduces notions of dissipa- 
tivity and stability for this class of quantum systems. In this 
paper, we build on the results of [16] to obtain robust stability 
results for uncertain quantum systems in which the quantum 
system Hamiltonian is decomposed as H = H\ + Hi 
where H\ is a known nominal Hamiltonian and H2 is a 
perturbation Hamiltonian, which is contained in a specified 
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set of Hamiltonians W. The set of perturbation Hamiltonians 
W corresponds to the set of exosy stems considered in [16]. 

For this general class of uncertain quantum systems, a 
number of stability results are obtained. The paper then 
considers the case in which the nominal Hamiltonian Hi is 
a quadratic function of annihilation and creation operators 
and the coupling operator vector is a linear function of 
annihilation and creation operators. This case corresponds 
to a nominal linear quantum system; e.g., see [3], [4], [8], 
[10], [15]. In this special case, robust stability results are 
obtained in terms of a frequency domain condition. 

The remainder of the paper proceeds as follows. In Section 
ILTl we define the general class of uncertain quantum systems 
under consideration. In Section [Hi] we consider a special 
class of quadratic perturbation Hamiltonians and obtain a 
robust stability result for this case. In Section [TV] we consider 
a general class of non-quadratic perturbation Hamiltonians. 
In Section [Vj we specialize to the case of a linear nominal 
quantum systems and obtain a number of robust stability 
results for this case in which stability conditions are given 
in terms of a strict bounded real condition. In Section |VT] 
we present an illustrative example and in Section IVIII we 
present some conclusions. 

II. Quantum Systems 

We consider open quantum systems defined by parameters 
(S,L,H) where H = R x + H 2 ; e.g., see [12], [16]. The 
corresponding generator for this quantum system is given by 



g(x) = -i[x,H] + c(x) 



(i) 



where C{X) = \tf[X,L] + \[L\X)L. Here, [X,H] = 
XH — HX denotes the commutator between two operators 
and the notation ' denotes the adjoint transpose of a vector 
of operators. Also, Hi is a self-adjoint operator on the under- 
lying Hilbert space referred to as the nominal Hamiltonian 
and H2 is a self-adjoint operator on the underlying Hilbert 
space referred to as the perturbation Hamiltonian. The triple 
(S,L,H), along with the corresponding generators define 
the Heisenberg evolution X(t) of an operator X according 
to a quantum stochastic differential equation; e.g., see [16]. 
The problem under consideration involves establishing 
robust stability properties for an uncertain open quantum 
system for the case in which the perturbation Hamiltonian is 
contained in a given set W. Using the notation of [16], the set 
W defines a set of exosystems. This situation is illustrated 
in the block diagram shown in Figure Q] The main robust 
stability results presented in this paper will build on the 
following result from [16]. 
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Fig. 1. Block diagram representation of an open quantum system interacting 
with an exosystem. 



Lemma 1 (See Lemma 3.4 of [16].): Consider an open 
quantum system defined by (S, L, H) and suppose there ex- 
ists a non-negative self-adjoint operator V on the underlying 
Hilbert space such that 



g(V) + cV < A 



(2) 



where c > and A are real numbers. Then for any plant 
state, we have 

(V{t))<e- ct (V) + - 7 Vt>0. 
c 

Here V(t) denotes the Heisenberg evolution of the operator 
V and (■) denotes quantum expectation; e.g., see [16]. 

A. Commutator Decomposition 

Given set of non-negative self-adjoint operators V and real 
parameters 7 > 0, S > 0, we now define a particular set of 
perturbation Hamiltonians VVi. This set W\ is defined in 
terms of the commutator decomposition 



[V,H 2 ] = \y,z ] ]w-vJ\z,V] 



(3) 



for FeP where w and z are given vectors of operators. 
Here, the notation [z, V] for a vector of operators z and 
a scalar operator V denotes the corresponding vector of 
commutators. Also, this set will be defined in terms of the 
sector bound condition: 

1 



w'w< 



,t. 



S. 



T 



Indeed, we define 



Wi 



H 2 :3i 



such that (|U is satisfied 



and ® is satisfied W £ V 



(4) 



(5) 



Using this definition, we obtain the following theorem. 

Theorem 1: Consider a set of non-negative self-adjoint 
operators V and an open quantum system (S, L, H) where 
H = Hi + H 2 and H 2 £ Wi defined in ©. If there exists 
a V E V and real constants c > 0, A > such that 

- i[V, H^ + C(V) + [V, zJ] [z, V] + \z^z + cV < A, (6) 



(V(t))<e- ct (V) + ^^, Vt>0. 



Proof: Let V E V be given and consider Q(V) defined in 
O- Then 

G{v) = -i^iJii + ^yj-j^zV + ^k^ 

(7) 

using (|3). Now since V is self-adjoint [V, z^jt = [z,V]. 
Therefore, 



< (\y,zi] 



IW' 



([V,J] 



iw' 



= [V, z f ] [z, V] + i[V, z^w - iuS [z, V] + w^w. 
Substituting this into (0, it follows that 

G(V) < -J[V,ffi]+£(y)-(-[y,z t ][z,V] + ^2; t z + <J 

T 

(8) 

using (0]i. Hence, (|6]l implies (O holds with A = A + S. 
Therefore, the result follows from Lemma [TJ □ 

B. Alternative Commutator Decomposition 

We now consider an alternative set of perturbation Hamil- 
tonians W2. For a given set of non-negative self-adjoint 
operators V and real parameters 7 > 0, Si > 0, S2 > 0, this 
set W2 is defined in terms of the commutator decomposition 

[V,H 2 ] = [V,z]wl-wi[z*,V] 

+^[z,\y,z]]wi-±w2[z,\y,z]]* (9) 

for V E V where w\, W2 and z are given scalar operators. 
Here, the notation * denotes the adjoint of an operator. Also, 
the set W2 will be defined in terms of the sector bound 
condition 



1 



and the condition 



wiwl < —^zz* + Si 
T 



W 2 Wn < S 2 - 



(10) 

(11) 



Then, we define 

H 2 : 3wi , W2, z such that ( fTUb and dTTb 
are satisfied and (O is satisfied W E V 

•(12) 
Using this definition, we obtain the following theorem. 

Theorem 2: Consider a set of non-negative self-adjoint 
operators V and an open quantum system (S, L, H) where 
H = Hi + H 2 and H 2 E W 2 defined in <H2}. If there 
exists a V E V and real constants c > 0, A > such that 
H = [z, [V, z]] is a constant and 

-i[V, H^ + C{V) + [V, z][z*, V] + \zz* + cV < A. 

7 

(13) 



Then 



(V(t))<e- ct (V) + 



X + Si +nn*/4 + S 2 



Vt>0. 



Proof: Let V <G V be given and consider Q(V) defined in 
©. Then 

g(v) = -ttv.jrii + rcvj-tty^K + itui^.v] 

—i[i,w* 2 + iw2jJ,* (14) 



using ©. Now [V, z]* 
self-adjoint. Therefore, 



z*V -Vz* 



[z*,V] since V is 



o < (\y,z]-iwi)(\y,z]-iwi)' 
= \y,z][z*,v] + i\y,z]wi 

—iw\\z*, V] + wiwl 
and hence 

-i[V,z]wl +iwi[z*,V] < [V,z][z*,V] +wiw*. 
Also, 



(15) 







1 



-Li - IW 2 



1 



-11 - IW-2 



rMM 



-w 2 fJ- 



-LIW2 + W2W2 



and hence 



-W2H- 



-fm>2 < -/i/i* +W2?«2- 



2 r 2 r z - 4' 
Substituting (Q3]i and ([T6]l into ([T4]i, it follows that 

G(V) < -i[V,H 1 ]+L{V) + [V,z\[z*,V] 

H — -zz* + Si + nil* /4 + <5 2 

7 2 



(16) 



(17) 



using ( fTUb and ( fTTT i. Then it follows from ( fT3T > that 

0(V) + cV < A + Si + mi* /A + <5 2 . 
Then the result of the theorem follows from Lemma [TJ □ 

III. Quadratic Perturbations of the Hamiltonian 

In this section, we consider a set W3 of quadratic pertur- 
bation Hamiltonians of the following form 



H2 = \[C C T ]A 






where A <G 



is a Hermitian matrix of the form 



A 



Ai A 2 
Af Af 



(19) 



and Ai = A\, A 2 = Aj. Also, ( = E x a + E 2 a*. Here 
a is a vector of annihilation operators on the underlying 
Hilbert space and a# is the corresponding vector of creation 
operators. Also, in the case of matrices, the notation t refers 
to the complex conjugate transpose of a matrix. In the case 
vectors of operators, the notation # refers to the vector of 
adjoint operators and in the case of complex matrices, this 
notation refers to the complex conjugate matrix. 



The annihilation and creation operators are assumed to 
satisfy the canonical commutation relations: 
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where J 



I 

-/ 



= J 



e.g., see [9], [13], [15]. 



(20) 



The matrix A is subject to the norm bound 

2 



IAIK 



7 



(21) 



where || ■ |j denotes the matrix induced norm (maximum 
singular value). Then we define 



m 



H 2 of the form ( flST l such that 
conditions (fT9land (fJJJi are satisfied 



(22) 



Using this definition, we obtain the following lemma. 
Lemma 2: For any set of self-adjoint operators V, 



W3CW1. 
Proof: Given any _ff 2 € W3, let 



Ai 
A # 



A 2 
Af 






AiC + A 2 C # 
AfC + AfC # 



and 



c 



E\ E<> 



Hence, 



H 2 



w* z 



E* 



Ef 



a 



= E 



a 

a* 



(23) 



]E*AE 



a 



Then, for any V £V 

1 



[V, z V 



yC f AiC + ^C f A 2 c # 

-V( T A* C + V( T A* (* 
CVA 1 c + C t ^A 2 c # 
+C T FA 2 # c + C T VA* (* 



(18) Also, 



•V\ = o 



C f AiCV + C T Af C^ 
-CtA 2 c # y + C T Afc # ^ 
1/ C f ^A!C + C T ^A#c 

2 V +Ct^A 2 C # + C T ^Af C # 
Hence, 

[V,z^]w-w^[z,V] 

= \ [V0 AiC + ^C f A 2 C # + ^C T A 2 # C + VC T A*Q* 

1 / C f AiC^ + C T Afc^ 



2 V +C f A 2 C # ^ + C T AfC # y 
Vff 2 - H 2 V = [V, #2] 



and thus (01 is satisfied. Also, condition (|2TT > implies 



t[(t C T ]AA 






< 



i[c f c T ] 



7 






which implies (@|i for any 5 > 0. Hence, H2 <G Wi. Since, 
Hi € W3 was arbitrary, we must have VV3 C W\. □ 

IV. Non-quadratic Perturbation Hamiltonians 

In this section, we define a set of non-quadratic pertur- 
bation Hamiltonians denoted W4. For a given set of non- 
negative self-adjoint operators V and real parameters 7 > 0, 
Si > 0, §2 > 0, the set W4 is defined in terms of the 
following power series (which is assumed to converge in 
the sense of the induced operator norm on the underlying 
Hilbert space) 



00 00 



00 00 



h 2 = /(c, = 2^2^ Skiver = i_i^ s U H U . 

fc=o £=0 fc=o e=o 

(24) 

Here S^e = Sj! k , H M = ( k ((*) e , and ( is a scalar 

operator on the underlying Hilbert space. It follows from 

this definition that 

OO OO DO OO 



fc=0 1=0 



e=o k=o 



and thus H2 is a self-adjoint operator. Note that it follows 
from the use of Wick ordering that the form (f24) is the most 
general form for a perturbation Hamiltonian defined in terms 
of a single scalar operator £. 
Also, we let 



00 00 



fc=i i=o 



fe-i 



(CY, 



00 00 



/" (c, n = EE fc ( fc - i)^c fc - 2 (o' 
fc=i £=0 

and consider the sector bound condition 

/'(c,c)7'(c,c)<^cr+^ 
7 

and the condition 

f"(cx*yf"(c,n<52. 

Then we define the set W4 as follows: 

H 2 of the form (|24j such that 



W 4 = 



conditions (l27l i and (f28t are satisfied 



(25) 



(26) 



(27) 



(28) 



(29) 



Note that the condition d28t effectively amounts to a global 
Lipschitz condition on the quantum nonlinearity. 

In this section, the set of non-negative self-adjoint opera- 
tors V will be assumed to satisfy the following assumption: 

Assumption 1: Given any V E V, the quantity 

A* = [C, [V, C]] = C[V, C] - [V, C]C 

is a constant. 



Lemma 3: Suppose the set of self-adjoint operators V 
satisfies Assumption [TJ Then 

W 4 c W 2 . 
Proof: First, we note that given any V E V and k > 1, 

VC = [VX]+CV; 



v( k = ^c"" 1 [^C]C fc "" + C fc ^- (30) 

n = l 

Also for any n > 1, 

C[V,C] = [v,C]C + w 



C-M^C] - [V.Cir-' + Cn-l^-V (3D 
Therefore using (l30l and (l3TT l. it follows that 
fe 

v^c fc = 53 [v, c]c n_1 c fc_ " + (n - i)c"- 2 c fc ~'v 

+C fc v 



fc[v,C]C 



fe-i + fc ( fc ^ c^/i + C^ 



which holds for any k > 0. Similarly 

(o fe ^ = fc(o fc ~ 1 [e*,yi + fc(fc 2 ~ - M*(O fc ~ 2 

+v(C) k - 

Now given any i?2 S W4, fc > 0, ^ > 0, we have 

[v,h m ] = kivxK'-'iCY + ^Y^K k - 2 (CY 
+C k v(cY 

-k( e (cY^[C,v} - fe(fc ~ 1 V c*(C*) fc ~ 2 
-CV(C) fc 
*[v,c]C*~ :l (C*) < -*C'(C*)*~ 1 [C*,v] 

) MC *- 2(0 / 



fc(*-l)..fc_2^ 



_fc(fc_ i )^ c£(c)fe _ 2 _ (32) 

Therefore, 

00 OO 

fe=o e=o 

= [v,C]f'{C,C)-f'te,C)*[C*,v] 

+^/"(c,c)-^*/"(c,cr. (33) 

Now letting 

z = C, tei=/ , (C,CT. andt« 2 =/"(C l CT. (34) 



it follows that condition (|9]l is satisfied. Furthermore, con- 
ditions ( [Tol l. ( fTTb follow from conditions ( |27| |. (l28l respec- 
tively. Hence, i?2 € W2. Since, 7^2 G W4 was arbitrary, we 
must have W 4 C yV 2 - D 

V. The Linear Case 

We now consider the case in which the nominal quantum 
system corresponds to a linear quantum system; e.g., see [3], 
[4], [8], [10], [15]. In this case, we assume that Hi is of the 
form 



1 



* = £[ 



,t 



M 



a 

a* 



(35) 



where M G 



is a Hermitian matrix of the form 



M 



Mi M 



M* 



M 



±2 
# 



ut 



and Mi = M-/ , M2 = M2 . In addition, we assume L is of 
the form 

L = \ JVk iV 2 



where JVi G 

L 

L# 



N 



and A^ 2 € 



a 



« 



Also, we write 



(36) 



N x N 2 



T# 



NT Nl 



it 



a 

n* 



In addition we assume that V is of the form 

V = [ at a T ] P 



a 

n* 



(37) 



where P G 
of the form 



•>2nx2n 



is a positive-definite Hermitian matrix 



P = 



P Pi 



P* 



P* 



(38) 



Hence, we consider the set of non-negative self-adjoint 
operators V\ defined as 

V of the form ([37) such that P > is a 



7>1 = 



Hermitian matrix of the form 



(39) 

In the linear case, we will consider a specific notion of 
robust mean square stability. 

Definition 1: An uncertain open quantum system defined 
by (S, L, H) where H = Hi + H 2 with Hi of the form 
(123, H 2 £ W, and i of the form (l36l l is said to be robustly 
mean square stable if for any H2 <G W, there exist constants 
c\ > 0, c 2 > and C3 > such that 
t r 



a(t) 

a*{t) 



< 



C\& 



-at 



a(t) 

a*{t) 

a 

a* 



a 

n* 



c 3 Vi > 0. (40) 



Here 



a(t) 
a*(t) 



denotes the Heisenberg evolution of the 



vector of operators » ; e.g., see [16]. 

In order to address the issue of robust mean square 
stability for the uncertain linear quantum systems under 
consideration, we first require some algebraic identities. 



Lemma 4: Given V € V\, H\ defined as in d35l ) and L 
defined as in ( l36l l. then 



[V,Hi 



[at a T ] P 



a 

n* 



,i[at o?]M 



a 



(I 



[PJM - MJP] 



a 



Also, 
C(V) = 



\tf[V,L] + \p,V]L 



Tr ( P JN f 



I 




1 

Furthermore, 



a* 



it 



NJ 

(nUnjp + pjn^jn) 



a 

a* 



a 



[ 



]P 



a 

a* 



2JP 



a 

,^- 



Proof: The proof of these identities follows via straightfor- 
ward but tedious calculations using (f2Qb . □ 

A. Quadratic Hamiltonian Perturbations 

We now specialize the results of Section HHIto the case of a 
linear nominal system in order to obtain concrete conditions 
for robust mean square stability. In this case, we use the 
relationship d23l : 






Ei E2 



r# 



E? El 



ii 



a 



E 



a 

a* 



(41) 

to show that the following following strict bounded real 
conditions provides a sufficient condition for robust mean 
square stability when H 2 G W3: 
1) The matrix 



F 



UM - -JN^JN is Hurwitz; (42) 
2 



2) 



E(sl-F) l D 



< 



1 



(43) 



where D = iJE^. 

This leads to the following theorem. 

Theorem 3: Consider an uncertain open quantum system 
defined by (S,L,H) such that H = Hi + H 2 where Hi 
is of the form d35j, L is of the form (O and H 2 G W 3 . 
Furthermore, assume that the strict bounded real conditions 
d42l . (l43l are satisfied. Then the uncertain quantum system 
is robustly mean square stable. 

Proof: If the conditions of the theorem are satisfied, then it 
follows from the strict bounded real lemma that the matrix 
inequality 



Ftp + PF + APJE^EJP 



E~<E 



<0. 



(44) 



will have a solution P > of the form d38k e.g., see [2], 
[4]. This matrix P defines a corresponding operator V G V\ 
as in ( |37V Now it follows from Lemma [4] that 



provides a sufficient condition for robust mean square sta- 
bility when H2 € W4: 

1) The matrix 



[z,V] 



2EJP 



a 



(45) 



F 



-iJM - -JN^JN is Hurwitz; (48) 
2 



where z is defined as in d23l i and d4TV Hence, 



2) 



[V,zt][^V]=4 






PJE^EJP 



a 



Also, 



E*Y,{sI-F) 1 D 
where D = JY,E T and S 



< 



(49) 



z^z = 



a 



B^ 



Hence using Lemma |4j we obtain 



a 



z* z 



-i[V, H^ + C{V) + [V, z ] ] [z, V] + - T + cV 



T 



7# 



F'P + PF 
APJE^EJP + 2lf 



a 

nit 



Tr PJN^ 



I 




NJ 



(46) 



where F = —iJM — ^JN^JN. Therefore, it follows from 
(l44t that there exists a constant c > such that condition 
(O will be satisfied with 



A = Tr PJN 



r1 



NJ >0. 



/ 


Hence choosing <5 = 0, it follows from Lemma [2] Theorem 
U and P > that 






a(t) 
a#(t) 



< 



e 



e" 



TO) 

(V) 



A 



0(0) I 1 1" a(0) 

a#(0) a # (0) 



= [-P] 



A 



CA-mi-n [r J 



v*>o. 



Hence, the condition d40b is satisfied with ci = 
c 2 = c> and c 3 = A , , > 0. 



i[P] 



_ JW;JP1 



7 
7 

This leads to the following theorem. 

Theorem 4: Consider an uncertain open quantum system 
defined by (S,L,H) such that H = Hi + H 2 where H x 
is of the form (l35l l, L is of the form d36t and 77 2 € W4. 
Furthermore, assume that the strict bounded real condition 
(|48T >. d49l is satisfied. Then the uncertain quantum system is 
robustly mean square stable. 

In order to prove this theorem, we require the following 
lemma. 

Lemma 5: Given any V € Pi, then 

H = [z, [z, V]} = [z*, [z*,V]}* = -ET.JP JE T . 

which is a constant. Hence, the set of operators V\ satisfies 
Assumption Q] 

Proof: The proof of this result follows via a straightforward 
but tedious calculation using < f2Qb - □ 

Proof of Theorem [?] If the conditions of the theorem are 
satisfied, then it follows from the strict bounded real lemma 
that the matrix inequality 

F ] P + PF + APJY,E T E*Y 1 JP+—Y 1 E T E*Y> < 0. (50) 

T 

will have a solution P > of the form d38l : e.g., see [2], 
[4]. This matrix P defines a corresponding operator V E Pi 
as in ( |37| |. 

It follows from ( f4Tb that we can write 



E*a* + E*a = [ E* E* 
E*T 



a 

n* 



AP\ 



> 0, Also, it follows from Lemma [4] that 



□ 



cA m ,„[P] 

B. Non-quadratic Hamiltonian Perturbations 

We now specialize the results of Section [IV] to the case of 
a linear nominal system. In this case, we define 

z = ( = E ia + E 2 a* 



[E, E 2 ] 



a 

n# 



E 



a 

n# 



(47) 



where z is assumed to be a scalar operator. Then, we show 
that the following following strict bounded real condition 



[**> 


v} = 


2E*Y,JP 


a 

. a * . 






lence, 

[V,z][z*,V} = 4 


a 

a* _ 


t 
PJY,E T E*T,JP 


a 

a* _ 


Uso, we can write 




zz — 


a 

. a * 


-it 

ZE T E* 
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Hence using Lemma [4] we obtain 
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where F = -iJM - \JNUN. 

From this, it follows using ( T50b that there exists a constant 
c > such that condition ( fT3l l is satisfied with 
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Hence, it follows from Lemma [5] Lemma [3] Theorem [2] and 
P > that 
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where A = A + 6\ + /U/i*/4 + fo- Hence, the condition d40b 



is satisfied with c\ 
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> 0, c 2 = c > and c 3 
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VI. Illustrative Example 
We consider an example of an open quantum system with 

S = I, #1=0, H 2 = -i (V) 2 - a 2 ) , L = V^a, 

which corresponds an optical parametric amplifier; see [17]. 
This defines a linear quantum system of the form considered 
in Theorem with M\ = 0, M 2 = 0, Ni = y^> ^2 = 

0, Pi = 1, P 2 = 0, Aj = 0, A 2 = i. Hence, M = 0, 
" Vk _ 
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E = T, and D = iJ. In this case, 
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which is Hurwitz, 
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VII. Conclusions 

In this paper, we have considered the problem of robust 
stability for uncertain quantum systems with either quadratic 
and non-quadratic perturbations to the system Hamiltonian. 
The final stability results obtained are expressed in terms 
of strict bounded real conditions. Future research will be 
directed towards analyzing the stability of specific nonlinear 
quantum systems using the given robust stability results 
for the case of non-quadratic perturbations to the system 
Hamiltonian. 
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Hence, we can choose 7 = 1 to ensure that (t2TT > is satisfied 
and H 2 G W3. Also, 
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Hence, it follows from Theorem [3] that this system will be 
mean square stable if - < |; i.e., n > 4. 



